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THE CRACK-INCLUSION INTERACTION PROBLEM* 

by 

Liu Xue-Hui and F. Erdogan 
Lehigh University, Bethlehem, PA 


ABSTRACT 


In this study the general plane elastostatic problem of interaction 
between a crack and an inclusion is considered. The Green's functions for 
a pair of dislocations and a pair of concentrated body forces are used to 
generate the crack and the inclusion. The integral equations of the problem 
is obtained for a line crack and an elastic line inclusion having an arbi- 
trary relative orientation and size. The nature of stress singularity 
around the end points of rigid and elastic inclusions is described. A ques- 
tion of specific interest which is studied is the nature of stress singu- 
larity around the point of intersection of the crack and the inclusion. 

Three special cases of this intersection problem which have been studied 
are a crack and an inclusion which are col linear and have a common end point, 
a crack perpendicular to an inclusion with a common end point (the L con- 
figuration), and a crack perpendicular to an inclusion terminating at its 
midpoint (the T configuration). The problem is solved for an arbitrary 
uniform stress state away from the crack-inclusion region. First, the non- 
intersecting crack-inclusion problem is considered for various relative size, 
orientation, and stiffness parameters and the stress intensity factors at 
the ends of the inclusion and the crack are calculated. Then for the crack- 
inclusion intersection case special stress intensity factors are defined 
and arc calculated again for various values of the parameters defining the 
relative size and orientation of the crack and the inclusion and the stiff- 
ness of the inclusion. 


T*T 

' ^ This work was supported by NASA-Langley under the Grant NGR 39 007 Qll 
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1. Introduction 


In studying the fracture of multi -phase materials, structures composed 
of bonded dissimilar solids, and welded joints it is necessary to take into 
account the effect of the imperfections in the medium. Generally such imper- 
fections are in the form of either geometric discontinuities or material 
inhomogeneities. For example, in welded joints various shapes of voids, 
cracks, notches and regions of lack of fusion may be mentioned as examples 
for the former and variety of inclusions for the latter. From the viewpoint 
of fracture mechanics two important classes of imperfections are the planar 
flaws which may be idealized as cracks and relatively thin inhomogeneities 
which may be idealized as flat inclusions with "sharp" boundaries. In both * 
cases the edges of the defects are lines of stress singularity and, conse- 
quently, regions of potential crack initiation and propagation. 

The technical literature on cracks, voids and inclusions which exist 
in the material separately is quite extensive. However, the problems con- 
cerning the interaction of cracks, voids and inclusions do not seem to be 
as widely studied (see, for example, [1] for the results of crack-circular 
inclusion or void interaction problem and for some references). In this 
paper the relatively simple problem of an elastic plane containing a crack 
and an arbitrarily oriented flat elastic inclusion is considered. Of special 
interest is the examination of the asymptotic stress field in the neighborhood 
of inclusion ends and the problems of intersecting cracks and inclusions. 

The basic dislocation and concentrated force solutions are used to formulate 
the problem [2]. Hence, the formulation can easily be extended to study prob- 
lems involving multiple cracks and inclusions. 

2. Integral Equations of the Problem 

The geometry of the crack-inclusion interaction problem under considera- 
tion is shown in Figure 1. It is assumed that the medium is under a state of 
plane strain or generalized plane stress and the ih-plane dimensions of the 
medium are large compared to the lengths of and the distance between the 
crack and the inclusion so that the effect of the remote boundaries on the 
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perturbed stress state may be neglected. Thus, the Green's functiohs for* 
the concentrated forces and dislocations in ah infinite plane may be u^ed tc^ 
formulate the problem. It is fuHher assumed that the inclusion is suffi- 
ciently "thin" so that its bending stiffness may also be neglected. 

Referring to Figure 1 we consider* the stresses and displaceihents due to 
a pair of edge dislocations on the x axis, a pair of concentrated forces bn 
the line 9=®c0nstant and the applied loads acting on the medium away from the 
crack-inclusion region. Let the subscripts d, p and a designate these thiree 
stress and deformation states, i.e., let and (iij) » (xiir) 

or (i,j) = (r,e), be the stress components due to dislocations, concentrated 
forces, and applied loads, respectively, the total stress state in the elas- 
tic plane may, therefore^ be expressed as 

cr^j(x,y) = oj^j(x,y) + Opjj(x,y.) + ag^.j(x,y)i (i,j = x,y) (1) 


Let us now assume that the dislocations are distributed along a-<x<bj 
y=0 forming a crack. If g(x) and h(x) refer to the dislocation densities 
defined by 


[Uy(x,+0) - Uy(x,-0)] = g(x) , a < x < b , 
^ [u^(x,+0) - u^(x,-0)] = h(x) , a < X < b , 


(2a, b) 


the corresponding stress coirponents at a poi nt (x,y) in the plane may be 
expressed as 




I CGjjj((x,y,t)g(t) + Hjjj,(xiy,t)h(t)]dt , 
a 




a 

rb 


[Gyy(x,y,t)g(t) + Hyy(x,y,t)h(t)]dt , 


[G^(x,y,t)g(t) + H^(x,y,t)h(t)]dt , 


(Sa-c) 

I 

i 


where 
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In (4) y and k are the elastic constants of the medium, y the shear modulus, 
K = 3-4v for plane strain and < = (3-v)/(l+v) for plane stress v being the 
Poisson's ratio. 

Similarly, from the concentrated force solution as given, for example, 
in [2] the stress components due to a pair of forces and P^ 

acting at the point (x^, y^) may be written as 


(A^+A2)Px + (B^+B2)P^ 
lCx-x„)i^ + (y-yj^^ 


Syy(x.y.Xo.yo) = 


(A^-A2)Pj^ + (B^-B2)P 


(5a-c) 


S^(x.y.x„.y„) = 25 ^ 


¥x ^ 63P. 


A, = -2(x-x„)[(x-x„)A + (y-y^)2] 

Ag = -<(x-Xq)[(x-Xp) 2 + (y-y^)2] - (x-Xq)[(x-Xq) 2 - (y-y^)2J + 2(y-yp)2(x-Xp) 
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Bi = •■2(y-yQ)C(x-x^)2 + (y-y^)2] 

B2 = +K(y-yQ)[(x-XQ)2+ (y-yQ)2]-(y-yQ)[(x-XQ)2-(y-y^)2]-2(x-x^j)2(y.y^) 
A3 = -K(y“yQ)C(x-XQ)2+(y-y^)2].(y>y^)[(x-XQ)2-(y-y^)2]^2(x-XQ)2(y-y^) 
83 = -<(x-XQ)C(x-XQ)2+(y-y^)2]+(x.XQ)C(x-XQ)2-(y-y^)2]-2(y-y^)2(x-j<^) 

(6a-^f) 

If the inclusion is located along the line c < > < d, e = constant, and 
if its bending stiffness is neglected, then the following conditions ane 
valid: 

Uy.(r*0+O) = u^(r,0-O), Ug(r,0+O) = Ug(r,0-O) , 

-PQ(r,0) = agg(r,0+O) - aQQ(r,0-O) = 0 , (7a-d) 

-Pp(i".e) = -p(r) = a^g(r,0+O) - a^Q(r,0-O) , (c<r<d). 

Thus, to forniulate the problem it is sufficient to consider only the radial 
component P^=p of the concentrated force. For Pg=0 and P^=p observing tha|: 

P^ = p COS0 . Py = p sin0 , (8a»b) 

and substituting x^ = r^cose, y^ = r^sina, by using the kernels given 
by (5) the stress components are found to be 

1 {A-| '+A2' )cos0+(B-j '+B^‘ )sin0 

^pxx^"^’-^^ " 27 t(ic+ 1) J L(^r^cos0)2+i:y-rQsin0)2i » 

c 

^ ^ 1 (A^ '-Ag' )cos0+(Bi '-B^' )sine 

apyy(x,y) - 2 tt(k+i ) J [(x-rQCOs0)2+(y-rQSin0)^J^ P^’^o^'^‘^0 * (Sa-c) 

^ ^ 1 A3' cose + B3'sine 

" 2ir(K+l) J Ux-l-^cos9)2+fy-r^sin9)2]2 P<''o>‘‘''o ’ 

Q 
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where the functions A^.', B/, (i=l,2,3) are obtained from (6) by substituting 
Xo = r^cose and = r^sine, e.g., 

A^'(x,y,rQ) = -2(x-rQC0se)[(x-rQC0se)2 + (y-rQSine)2] . (10) 

Since the stresses due to the applied loads are known, from (1), 

(3) and (9) it is seen that once the functions g(x), h(x) and p(r) are deter- 
mined the problem is solved. These unknown functions may be determined by 
expressing the stress boundary conditions on the crack surfaces and the 
displacement compatibility condition along the inclusion, namely 

ayy(x,0) = (Jjyy(x,0) + Opyy(x,0) + <Jayy(x,0) = 0 , (a<x<b), 

o^yi^yO) = cfjj^y(x,0) + «Jpxy(^»0) + ^ * (a<x<b), (lla-c) 

= edy.^(r,0) + ep^y.(r»e) + ^arr^*"’®^ e^-(r), (c<r<d) 

where £^(r) is the (longitudinal) strain in the inclusion. If, for example, 
the stress state away from the crack inclusion region is given by 
(i»j) = (x,y), then the applied quantities in (11) may be expressed as 


^ayy(^»°) = cr”y , cJg^y(x,0) = , 


earr^*"’®^ = ^ CO (cos2e - sin2e) 


8y '• XX 


I't'K 


+ - f?f * TT? C sin2e] 


(12a-c) 


We now note that if p(r) is the body force acting on the elastic medium 
then -p(r) would be the force acting on the inclusion distributed along its 
length. Thus, the strain in the inclusion may be obtained as 

d 


f 


031 


- 5 - 



where and ahe the elastic cohstarits, and is the cross-sectional 
area of the irtcliJsiori corhespdndiiici to lihit thickHdss bf the medium in z- 
direction. From the expressibri of given by the Hbbke's law 


c. = 3 -k 


rr 8y '^rr 1 +k 


^iaa ) » 


(14) 


frbm (9) and the cbrrespbhdiri^ sti^ss thahsl'bi’metibH it can be shown that 

d 




ku! 


PO^h) 


dr . . 


h.-r "0 


Similarly, fi"om (3), (4) ahd (14) Wb fihd 

b 

^drr^*^’®^ = W I + Hjr,t)h(t)]dt 


where 


® ^ Tr(l+K) ^cos 29 - 1^ sin2e)(t-r cose) x 

x[(t-rcose)2-r2sih2eJ 4- (sin^e- cos^e) x 

x(t-rcos0)[3r2siH20 + (t-rcose)2] 


(1^) 


(16) 


+ sin29 r 9irl0[r2siH28-(t-rcos0)2]} , 
H^(r,t) = ^ {(co 02 e- sin20)r sinei[r 2 sin 20 

+ 3(t-rCos0)2] + (sih20- cbs20)r sine x 
x[r 2 sirt 20 -(t-r 2 ccis 2 ei| 2 ] + ^ sin 2 e x 
x(t-hcose)[(t-i'CbSej2 - ^'2^in2ej , 


r 2 = (t-r cos 0)2 + r 2 si 


(17) 


(18) 

(19) 
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Finally, by substituting from (3), (4), (9), (12), (13), (15) and 
(16) into (11), the integral equations of the problem may be obtained as 
follows: 


1 f” q(t)dt , 1 (Ai'-A;')cose^(B^'-B 3 ')s 1 ne 

7T J t-x 47ry J [ { x-r COS 0 ) ( r s i n 0 ) ^ 2y ^yy ’ 


1 [ h(t)dt . 1 f (*3'c°se-i-B3'sin6)p(r^) 
tr J t-x 4 irp J [(x-rQCOS9)^+(r^sine)^]^ ^^'o 


(a<x<b) , 

1+K 


2u "xy ’ 


_o 

IT 


r C /•“ , 

G_(r,t)g(t)dt + — H (r,t)h(t)dt + - 

J S IT j S 7 T 


(a<x<b) , 

— dr 

r^-r 0 


~T I [(cos20- sin20)a 


00 

XX 


+ (sin20- cos20)0”y + ^ sin2o], (c<r<d) , 


where 


c„ = 


41+k)2 . “<’%> 

4 ic ’ ■ AjUjO+k) 


From the definition of g and h given by (2) it follows that 
b b 


g(t)dt = 0 , 


h(t)dt = 0 . 


Also, the static equilibrium of the inclusion requires that 


d 

* 

p(r)dr = 0 . 

« 

c 


(2da-c) 


(21a, b) 


(22a,b) 


( 23 ) 
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Thus, the system of singular integral equations must be solved under the 
conditions (22) and (23). From the function-theoretic examination of the 
integral equations (20) it can be shown that the unknown functions g, h and 
p are of the following form [2]: 


g(t) = 


Fi(t) 


(b-t)i(t-aji 


h(t) = 


F,(t) 


(b-t)^(t-a)^ 


P(r) = 


F3(r) 


(d-r)^(r-c)^ 


(24a-c) 

where F^, and F^ are bounded functions. The solution of (20) subject to 

(22) and (23) may easily be obtained by using the numerical method described 
in [3], 


3. Stress Singularities 

After solving (20) the Modes I and II stress intensity factors and 
k 2 at the crack tips x=a and x=b, y=0 which are defined by 


k,(a) = lim »^2(a-x) a (x,0) , k,(b) = lim /2(x-b) a,„Tx,0) , 

' v.vo jj I .. yy 




x-*b 


kp(a) = lim y2(a-x) 0 ^„(x,O) , 
x->a ^ 


may be obtained as follows: 


kp(b) = lim v^2(x-b) ct„,,(x,0) , 
^ x->b 


(25a-d) 


ki(a) = ^ lim i/ 2 (x-a) g(x) , k,(b) = - lim i/ 2 (b-x) g(x) , 

x->a x-it) 

•<2(3) - ^ li'T’ *^^(x-a) h(x) , k2(b) = - lim / 2 (b-x) h(x) . 

x->a x->b 


(26a-d) 


The constants k-j and k 2 related to the asymptotic stress fields near 
the crack tips through the well-known expressions (see, for example, [4] 
and [5]). However, not so well-known is the asymptotic behavior of the 
stress fields near the inclusions having sharp edges. From (24c) and (7d) 
it is seen that the shear stress has a square-root singularity at the 
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tip of the inclusion. However, if one is interested in crack initiation 
around such singular points, one needs to know the direction and the magni- 
tude of the maximum local cleavage stress. This, in turn, requires the 
investigation of the complete asymptotic stress field near the singular 
points. By using the basic form of the solution of the related density 
functions given by (24) and going back to the original stress expressions, 
the asymptotic stress fields may be developed by following the general 
techniques described in, for example, [6] or [7]. 

In an elastic medium containing an elastic line inclusion under plane 
strain or generalized plane stress conditions, the asymptotic analysis gives 
the near tip stress field as follows £7]^*^ 



3+k 

K-1 





K+l 

K-1 



£ 

2 ’ 


(27a-c) 


where x,y and r,e are the standard rectangular and polar coordinates, the 
origin of coordinate axes is at the inclusion tip and the inclusion lies 
along the negative x axis or along e=,r, r>0. Equations (27) suggest that 
similar to crack problems one may define a (Mode I) "stress intensity factor" 
in terms of the (tensile) cleavage stress as follows: 


k^ = lim »^ayJr,0) . 


(28) 


From (7) by observing that (at the right end of the inclusion) 

axy(r.+TT) - a^y(r,-ir) = -p(r) , (29) 

7*y 

Note the misprints in (4.6) of £7]. 
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(30) 


in terms of the function p(x) may be expressed as 


k, = 


•lim 

r->0 




It should be noted that in the case of flexible elastic line inclusions 
there is no antisymmetric singular stress field. For example, in a plane 
under pure shear (<y^y) parallel to the inclusion, the perturbed stress field 
is zero. Physically this of course follows from the fact that the normal 
strain parallel to the plane of shear is zero. 

Similarly, for a rigid line inclusion (i.e., for an inclusion having 
infinite bending as well as tensile stiffness) it can be shown that for 
small values of r the asymptotic stress field Is given by 


^yy(r.e) = 


1 


(ki cos f + ^ k, sin f) 


cTx„(r,e) s _L (_ ItK k 

XX ^ ' K-1 1 


cos ^ + 


0 . 3-ic 


K-1 


k 2 sin 


2 ^ 


(31a-c) 


Again, the stress intensity factors k^ and k 2 are defined in terms of the 
tensile and shear cleavage stresses at e=0 plane as follows: 


k^ = lim y(r,0) , k 2 = lim v^a^^(r,0) . (32aib) 

As in the crack problems, the anti plane shear component of the asymp- 
totic stress field around flat elastic and rigid inclusions is uncoupled. 
Defining a Mode III stress intensity factor by 


k 3 = lim *^ax-(r,0) , 
r-K) 

the asymptotic stress field may be expressed as 


(33) 
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(34a,b) 


•'O n 

o^y{r,Q) = — ^ cos f- , 

~ A 

a (r,e) = — sin j , 

where again the inclusion lies along e=ir plane^*\ 

4. Crack-Inclusion Intersection 

Analytically as well as from a practical viewpoint intersection of 
cracks and inclusions presents some interesting problems. In these problems 
the point of intersection is a point of irregular singularity with a power 
other than 1/2. Even though the general intersection problems for an arbi- 
trary value of 9 may be treated in a relatively straightforward manner, in 
this paper only some special cases will be considered. 

4.1 The case ofe=|-,a-0, c = 0 


In this case the system of singular integral equations (20) becomes 


i 

TT 


1 

IT 


f aftl 1 r 

^ + - J - ■(■j+^ 2 ) 2 ]p(t)dt = f^(x) , (0<x<b) , 


0 

b 


0 

d 


^ I Ccxi+t2)2 ” ; 5 ^ 2 i^]p(t)dt = f 2 (x) , (0<x<b) , 


TITT " 

'Note that in this case if the remote stress is decomposed into a^z and 
0 “^, the perturbed stress field due to 0 “^ would be zero. For the cleavage 
plane 0 the shear cleavage stress may be written as 0 - {r, 8 )= 0 „ sine -0 c6se= 

vQ XZ 

-(k 2 /v?r)sin(e/ 2 ) , = 0+ir/2, indicating that e='fTr/2 is the maximum cleavage 

planes. 
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0 0 



I 


p(t) 

t-r 


dt + 



ir 


d 

I H(t-r)p(t)dt = f^{r) , (0<r<d), 


where 






"3 


vU-l) 


4 < ’ ^5 4AgKyg 


(3Sa-c) 


(36) 


and f^, and are known input functions (see, for example, the right 
hand side of (20)). Note that aside from thp simple Cauchy kernels, (35) 
has kernels which become unbounded as the variables (t,x,r) approach the 
point of irregular singularity (x-0=t-r). Thus, defining the unknown func- 
tions by 


g(t) = 


Fi(t) 


h(t) = 


Fj(t) 


p(t) 


tv-t)®’ ’ 


0 <Re(a,e|^)<l, (k=l,2,3) , (37a-c) 

and by using the function-theoretic technique described in [ 3 ], the charac- 
teristic equations for 3^, Bg’ ^3 “ may be obtained as follows: 


COt:r3,^ = 0 , (k = 1,2,3) 

b-j cos^ira - (b 2 + 8 a-b 3 a^)cos 2 ^ 


-(b4-b5a+b3a2)sin2 ^ = Q , 


where 
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( 40 ) 


= 8k/(1+k) , 62 = 2{3+ic)(k:-1 )/(<+! ) , 

b3 = 8 /(k+ 1) , b^ = 2(3 -k) , bg = 16/(l+ic) . 

Note that the properties of the inclusion (as expressed by the constant Cg 
in (36)) enter the integral equations (35) only through a Fredholm kernel 
cind, therefore, have no influence on the singular behavior of the solution, 
and a is dependent on k or on the Poisson's ratio of the medium only. From 
(38) it is seen that the acceptable roots are = 0.5, (k = 1,2,3). The 
numerical examination of (39) indicates that in this special case of 9 = I" 
we have 0.5<a<l, meaning that the stress state at r=0=x has a stronger singu- 
larity than the conventional crack tip singularity of 1//F, This may be 
due to the fact that in this problem two singular stress fields are combined 
at r=0. Also, it turns out that for 0<v;^5 the characteristic equation (39) 
has two roots in 0<Re(a)<l and both are real. These roots are given in 
Table 1 for various values of the Poisson's ratio. 

Table 1. Powers of stress singularity a for a crack and 
an inclusion: a = 0, c = 0, 0 = ir/2 (Fig. 1). 



plane strain 



plane stress 



V 

“1 

“2 

“1 

“2 


0.63627093 

0 

0.63627093 

0 

0.1 

0.64489401 

0.09571474 

0.64408581 

0.08990596 

0.2 

0.65405762 

0.14825371 

0.65095281 

0.13249000 

0.3 

0.66352760 

0.18953334 

0.65695651 

0.16176440 

0.4 

0.67270080 

0.22567265 

0.66217253 

0.18404447 

0.5 

0.67996342 

0.26027940 

0.66666667 

0.20196313 


The stress intensity factors at the crack tip x=b, y=0 and at the end 
of the inclusion x=0, y=d may be obtained by using the relations (26) and ^ 
(30). At the singular point x=0, y=0 the following useful stress intensity 
factors are defined; 


- 1 3 -^ 



(41a, b) 


k,(0) = 11m x“ o (-0,0) , 
x-»-n yy 


l< 2 ( 0 ) = lim ^ x“ a (-0,0) , 

x-)-o ^ 

for the crack, and 

ki(0) = lim -5^y“ p(0,+0) 
y-vfO ^ 


for the inclusion. 


(42) 


4.2 The Special Case Of e = | , c =* -d, a = 0. 

In this case the problem is further simplified by assuming “symmetric" 
external loads (for example, a;^-0 in (20)). Thus, the plane of the crack 
is a plane of symmetry, h(x) = 0, and (20) would reduce to 


1 

7T 



2 C2tx2 

- J ^ * f^(x), (0<x<b) 

0 


9 


c "t c ty ^ ^ 

J J *^t2+y2 (t2+y2)2]g(t)0t +1 I ^ 


1 


+ GgH(t-y)ip(t)dt fg(y) , (0<y<d) , 


(43a, b) 


where, again the input 
given in (20) (with a" 
By defining 


functions f^ and fg are known and, for example, are 
= 0) and the constants c^ Cg’are defined by (36). 


g(t) = 


6^(t) 

t“(b-t)^i 


P(t) - 


G^it) 


(44) 


from (43) it may be shown that 
cote,^ = 0 , (k=l,2) , 
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(46) 


C0Sira-(C3 + ^ c^a)(c^ - j € 20 ) = 0 . 

From (45) it is seen that = 0.5. A close examination of (46) shows that 
it has only one root for which 0<Re(a)<l. Furthermore, this root turns out 
to be real and highly dependent on the Poisson's ratio (see Table 2). The 
characteristic equation (46) and the roots given in Table 2 are identical 
to those found in [8] where an infinitely long stringer in cracked plate 
was considered. 

Table 2. Power of stress singularity a at the crack-inclusion 
intersection for e=ir/2, c=-d, a=0 and for symmetric 
loading. 


V 

a 1 

plane strain 

plane stress 

0 

0 

0 

0.1 

0.10964561 

0.10263043 

0.2 

0.17432137 

0.15468088 

0.3 

0.22678790 

0.19132495 

0.4 

0.27392547 

1 0.21972274 

0.5 

0.31955800 

0.24288552 


In this problem, too, the stress intensity factors for the crack and 
the inclusion may be defined as in (41) and (42). 

4.3 The Special Case of 0=ir, a=0, c=0 

In this case the crack and the inclusion are on the x axis and occupy 
(y=0, 0<x<b) and (y=0, -d<x<0), respectively. Restricting our attention 
again to the symmetric loading for which h(x) = 0 and observing that for 
the variables along the inclusion r = -x, r^ = -t, p(r^) = -Pj((t), the 
integral equations of the problem may be expressed as 


/ 
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am 

t-x 



rr 


K-1 

4 m 


u 

I 


p„(t) 


f ^ dt + 

J t-x 
0 


1 f ^5 f 

- J dt - — j p^(t)dt = f3(x), (-d<x<0) (4^a.b) 


where the constants G3 and C5 are defined by (36) and the known functions 
f^ and f^ are given by the right hand sides of (20a) and (20c) (with c“ =0). 
If we now let 


g(t) = 


H^(t) 




Hg(t) 

(-t)“(t+d)®2 


0<Re(a,Bi,B2)<l . 


(48) 


from (47) the characteristic equations for a, 3i and &2 be Obtained as 
fol lows : 

cotirg,^ = 0 , (k = 1,2) , (49) 

COs2Tra = - (— )^ . (50) 

2/iT 

Equation (49) again gives = 0*5. From (50) it may easily be seen 

that a is complex and its value for which 0<Re(a)<l is found to be 

a = 1+ i(^^^) . (51) 

This value of a turns out to be identical to the power of singularity for a 
perfectly rough rigid stamp with a sharp corner pressed against an elastid 
half plane having < as an elastic constant [2] (e.g., < = 3-4v for the pl4ne 
strain case). At first this hesult may be somewhat unexpected* However^ 
upon closer examination of the problem first, from (47b) it may be seen 
that the elasticity of the inclusion (i.e., the term containing the constclnt 
Cg) has no effect on the nature of the stress singularity. Thus, if one 
assumes the inclusion to be inextensible, for the symmetric problem under 
consideration it can be shown that the conditions in the neighborhood of the 
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crack tip x=0, y=0, for example, for y<0, are identical to the conditions 
around the corner of the stamp in the elastic half plane occupying y<0. 

It, therefore, appears that for the elastic inclusion col linear with a 
crack, the stress state around the common end point would have the standard 
complex singularity found in the rigid stamp problem. 


5. The Results 

The crack-inclusion problem described in previous sections is solved 
for a uniform stress state a^.j, (i,j=x,y), away from the crack-inclusion = 
region. For simplicity the results are obtained by assuming one stress com- 
ponent (a“^ or a”y or a“^) to be nonzero at a time. The solution for a 
more general loading may then be obtained by superposition. Even though , 
the stress state everywhere in the plane can be calculated after solving the 
integral equations (e.g., (20)) and determining the density functions g, 
h, and p, only the stress intensity factors are given in this section. For 
nonintersecting cracks and inclusions the stress intensity factors defined 
by (26) and (28) are normalized as follows: 


^ k.(x.) 

x.=(a,b); 


J <j"vT152J)72 


a ' yy’ XX’ xy' 


(52) 


for the crack and 

■ . k,(r.) 




(rj = (c.d) , - (o" c" )) 


(53) 


for the inclusion. 

Referring to Figure 1, for c=a, d=b, and (b/a)=5 the effect of the angle 
9 on the stress intensity factors is shown in Table 3. These results are- 
given for two values of the stiffness parameter y defined by (21), namely 
Y=0 (the inextensible inclusion) and y=10. ' 

j 
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Table 3. Normal i zed stress intensity factors in a plane containing a 
crack and an inclusion subjected to a uniform stress state 
away from the crack*incl us ion region (c=a,d=b»a=b/5iFig. 1). 
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Table 3 - cont. 


m 

■■ 

9 

■ 



30® 

60® 

90® 

120® 

O 

O 

in 

O 
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00 

1 

■ 


(c) c 

00 , * 

’yy * ° ’ 

Q 

11 

O 

- r 

* ""xy ■ ^ 

) 




.1289 

.1428 

.0669 

.0028 

.0134 

.0223 

0.0000 


k^(a) 

1.0849 

1.0180 

.9054 

.9950 

1.0599 

1.0304 

1.0000 

0 

ki(b) 

.1641 

-.0754 

-.0670 

-.0021 

0.0231 

.01 36 

0.0000 


k^(b) 

1.4055 

.9685 

.9974 

.9995 

1.0005 

1.0005 

1.0000 


k;(c) 

-1.0246 

-1.6348 

-1.3085 

.0533 

1.3767 

1.3606 

0.0000 


ki(d) 

2.0539 

-1.3808 

-1.4661 

-.1076 

1.2735 

.3117 

.0000 


kl(a) 

WM 

.0198 

.0100 

.0010 

.0032 

.0043 



k^(a) 


.9967 

.9826 

.9992 

1.0108 

1.0054 


10 

k;(b) 


-.0140 

-.0121 

-.0003 

.0043 

.0025 

.0000 

k^(b) 

ng 

.9929 

.9994 

.9998 

.9999 

1.0000 

1.0000 


kUc) 


-.5492 

-.3731 

.0557 

.4513 

-.4316 

.0000 


k;(d) 

IB 

-.4179 

-.4533 

-.0342 

.3912 

-.4029 

.0000 
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Some sample results for an inclusion collinear with a crack (i.6, 
for 0=0) are given in Table 4. Note that for this configuration under the 


Table 4. Normalized stress intensity factors for an inclusion 
collinear with a crack. Relative dimensions: 6=0, 
d-c = b-a, c = b+s. Applied loads: c?., (i,j=x,y) 
(Fig. 1). 


00 

CJ • • 

k' 

s = (b-a)/100 

s = (b-a)/2 



Y = 0 

o 

II 

Y - 0 

Y = 10 


Bil 

-0.0202 

-0.0040 

-0.0019 


00 

rt 


-0.1338 

-0.0300 

-0.0027 


XX 

IH 

-1.0482 

-0.3296 

-1.0889 

BBI 


■B 

-1.0845 

-0.3345 

-1.0889 



ED 


1.0006 

1.0008 


00 

HD 


0.9987 

1.0011 


"yy 



-0.1571 

0.4559 

0.1397 



0.3841 

0.1273 

0.4590 

0.1413 


loads shown in the table, that is, for a“ and a”, because of symmetry 

yy XX 

the Mode II stress intensity factors k 2 (a) and k 2 (b) are zero. Also, for 
the shear loading it is found that k^Ca) = 1, k^Cb) = 1 and k-|(a) = 
k-|(b) = k-j(c) = k-j(d) = 0. This follows from the fact that in the cracked 
plane under pure shear a~ the strain component e^^(x,0) is zero and 
hence an inextensible inclusion on the x axis would have no effect on the 
stress distribution. 

Another special configuration is an inclusion parallel to the crack 
for which Table 5 shows some sample results. In the two special conv 
figurations considered in Tables 4 and 5 the effect of the crack-inclpsion 
interaction on the stress intensity factors does not seem to be very sig- 
ni f i cant . 

The results for an elastic medium for which xz plane is a plane of 
symmetry with respect to the crack-inclusion geometry as well as the 


■ /' 
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Table 5. Normalized stress intensity factors in a plane containing an 

inclusion parallel and equal in length to a crack, both symmetri- 
cally located with respect to the y axis. The crack is along the 
X axis and H is the distance between the crack and the inclu- 
sion in y direction (Fig. 1). 


oo 

"^ii 

k' 

H = 

b-a 

H = lO(b-a) 




Y = 10 

l■^ggQIIIIII 

O 

11 


k[(a)=l<2(b) 

-0.0182 

-0.0070 

-0.0007 

-0.0002 

00 

^x 

k'(a)=-k^(b) 


-0.0011 

0.0006 

0.0000 


ki(c)=kl(d) 


-1.0887 

-0.0683 

-0.0683 


k;(a)=ki(b) 

1.0063 



1.0001 

00 

"yy 

k^(a)=-k^(b) 

-0.0060 



0.0000 

ki(c)=ki(d) 

0.3917 

0.4387 

0.0411 

0.0276 


ki(a)=-ki(b) 

-0.0042 

0.0000 

-0.0002 


oo 

'^xy 

k2(a)=k2(b) 

0.9965 

1.0000 

0.9998 


k{(c) 

-0.1131 


-0.0123 

0.0004 


kl(d) 

0.1129 

-0.0052 

0.0123 

-0.0006 
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applied loads are given in Figures 2-12. In this example the crack is per- 
pendicular to the inclusion and the external load is a uniform tension par- 
allel or perpendicular to the crack and away from the crack-inclusion region 
(see the insert in the figures). The results shown in the figures are self- 
explanatory. However, the solution also has some unusual features among 
which, for example, one may mention the tendency of the crack tip stress 
intensity factors k'(a) and k'(b) to "peaking" as y decreases and as d/i 
increases (where 2d and 2i are the lengths of the inclusion and the crack, 
respectively and y * 0 corresponds to an inextensible inclusion). 

The results for the limiting case of the crack touching the inclusion 
are given in Figures 8-12. In this case at the singular point x=0, y^^O 
the stress intensity factor k-j(a) and the normalized stress intensity factor 
kj (a) are defined by 

(54) 

I 

(55) 

where the power of singularity a is given in Table 2. The results shown 
in Figures 8-12 are obtained for v = 0.3. 

The stress intensity factors for the other symmetric crack inclusion 
problem, namely for the problem in which y axis is the line of symmetry 
with regard to loading and geometry are given in Figures 13-28. In this 
problem a=-£, b=£, d>c>0 and the external load is either a” or a” (see the 

• yy XX 

insert in Figure 13). Note that the figures show the crack tip stress 

intensity factors at x=a=-£ and k^(b)=k^(a), k 2 (b)=-k 2 (a) . Generally the 
magnitude of k^(a) and k 2 (a) seem to increase with increasing length and 
Stiffness of the inclusion (i.e., with increasing (d-c)/2ji and decreasing 
Y = y(l+Kg)/Ag)i^(l+ic), where is the shear modulus of the inclusion). 

Also, as expected, k-|(c) and k-j(d) describing the intensity of the stress 
field at inclusion ends tend to increase as the stiffness of the inclusion 
increases. However, their dependence on the relative length parameters is 
somewhat more complicated (see, for example. Figure 16 for change in beha- 
vior of the variation of k^(d) at (d-c)/2ji=5). Figures 13-20 show the effett 


k-|(a) = lim ^/l x"^ (x,0) , (x < 0) , 

x->0“ 

k-j'(a) = k(a)/a?^/F , (i=(x,y); «,=b/2) 


-22r. 



of the inclusion length for constant crack length 2i and constant distance 
c (Figure 13). The effect of the distance c for constant inclusion and 
crack lengths is shown in Figures 13-28. 

The results of the nonsymmetri c problem showing the effect of the rela 
tive location of the inclusion are shown in Table 6. Referring to Figure 1 
in these calculations it is assumed that e = |- , d-c = 2t, c/2z = 0.1 and 
a/2z is variable. 

Finally, the stress intensity factors for the crack-inclusion inter- 
section problem considered in Section 4.1 are given in Figures 29-43. 

The normalized stress intensity factors shown in these figures are defined 
by (see (41), (52) and (53)) 



— lim /2(x-b) (T^^(x,0) , 

a”./Tx->b yy 

* J 


koD = — lim /2(x-b) a (x,0) , 

a . . /i x-+b y 

V w 

*^ 1 /\ “ ~Z P x“ cr (-0,0) , 

x->-o yy 

* \J 


1 


'2A 


x->-0 


lim ^ x“ _(-0,0) 


xy' 


_ 

'ID " k. 


knn = 7 ^ lim /2(y-d) a^^(0,y) 


0 y-Ki 


XX' 


(56) 

1 


In this case too, generally the magnitude of the stress intensity factors 
increases with increasing length and stiffness of the inclusion. However, 
since the crack and the inclusion are located in each other's "shadow", 
the relative dimensions seem to have considerable influence on the vari- 
ation as well as the magnitude of the stress intensity factors. 
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= 0. a/£> 0.5, v = 0.3. 
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Figure 6. Normalized stress intensity factor at the crack tip x=b: 0 




Normalized stress intensity factor at the inclusion end y=d; 
^ ^yy “ = 0* = 0-5. v = 0.3. 
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Normalized stress intensity factor at the crack tip x=a=0, 
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Figure 14. 
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Figure 15. 
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Figure 18. Mode II stress intensity factor at the crack tip x=a 
^ 0, a~y = cr~y = 0, V = 0.3, 9 = tt/2, c = 0.2il, 


n jo 




Figure 20. Stress Intensity factor at the inclusion end y = d; <r“ 0, 

CTyy - cr"y = 0, V = 0,3, 0 = ir/2, c = 0.2z, b = £ = -a^^ 












Figure 25. 
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Figure 26. Mode II stress intensity factor at the crack tip x = a = -li 
^ ^yy ~ ®xy * 0, v = 0.3, 0 * ir/2, d-c = 2t, b * A. 
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Figure 29. Normalized stress intensity factor for the inclusion-crack 
intersection problem for which 9 = ir/2, a = 0, b - 2a. c * 0 
d/2A and y variables, for / 0. = 0, 0 “^ = 0. ’ 
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Figure 30, Normalized stress intensity factor for the incl us ion -crack 

intersection problem for, which e = ir/2, a = 0, b = 2a, c = 0, 
d/2A and y variables, a” = 0, a" = a“ = 0. 
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Figure 32. Normalized stress intensity factor for the inclusion-crack 

intersection problem for, which 0 * ir/2, a = 0, b = 2z, c - 0, 
d/2ji and y variables, kgg . 5 ^ 0. 



Figure 33. Normalized stress intensity factor for the inclusion-crack 

intersection problem for, which e = ir/2, a = 0,b = 2jt,c = 0, 
d/2i and y variables, * 0. 





Figure 34. Normalized stress 
intersection probl 
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Figure 35. Normalized stress intensity factor for the inclusion-crack 

intersection problem for, which e = tt/ 2, a = 0, b = 2t, c = 0, 
d/2t and y variables, kgn , a“ t 0. 
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Figure 40. 
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